Internal waves and synchronized precession in a cold vapor 
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Exchange in a Boltzmann gas of bosons with several internal states leads to collective transport 
of internal polarization. The internal dynamics can be understood as Larmor precession in the 
presence of a torque induced by atoms on each other via exchange coupling. A generalized Bloch 
equation that includes interatomic exchange effects as well as orbital motion in the gas is derived 
and used to interpret recent experiment by Lewandowski et al. |l[ as an excitation of a collective 
wave of internal state polarization. It is shown that exchange leads to formation of domains in which 
precession frequencies are synchronized. 



Atomic gases in the cold collision regime characterized 
by de Broglie wavelength long compared to the range of 
interparticle potential represent an interesting quantum 
many-body system. Most surprisingly, spin waves in a 
cold spin polarized gas are collective excitations. This 
phenomenon was actively studied in the 80's, first pre- 
dicted by Bashkin and independently by Lhuillier and 
Laloe |?| , and confirmed by NMR experiments in spin po- 
larized H I by Johnson et al. Q, in 3 He by Nacher et al. 
|p|, and in dilute 3 He— 4 He mixtures by Gully and Mullin 
|3|. A detailed quantitative theory of the observed NMR 
spectra was given by Levy and Ruckenstein Q . Bigelow 
et al. Q demonstrated that collective spin waves are pre- 
served even in the Knudsen regime. The theory was fur- 
ther developed by Miyake et al. || and reviewed in jl0| . 

Exchange effects in gases are not limited to spin phe- 
nomena, since any pair of internal states can play a role 
similar to spin states in exchange collisions Apart 
from new energy scales arising due to internal states spec- 
trum, the main difference is in the anisotropic character 
of exchange. For generic internal states the Hamiltonian 
does not have spin-rotational symmetry, which opens a 
number of new interesting possibilities. Verhaar et al. 
|l3f demonstrated that interatomic exchange leads to en- 
hancement (by a factor of two) of the density shift of 
Rabi transition. Tiesinga et al. Jl4| and Kokkelmans et 
al. |l5| considered application of this effect in Cesium 
fountain clock. Similar exchange enhancement occurs in 
the optical spectrum density shift. The theory (^|| of 
collective spin-waves, can be straightforwardly adapted 
to describe optical excitations jj^ ]. 

New aspects of cold collision exchange arise in experi- 
ments on Bose-Einstein condensation (BEC) in trapped 
gases. The exchange part of the density shift is absent in 
BEC at T = and is reduced at < T < T BE c ■ 
Interestingly, in this case all modes involving coupling of 
internal states are split into doublets JlSj ] . 

In a recent experiment Rabi transition was studied 
in a cigar shaped sample of Rb vapor contained in Ioffc- 
Pritchard trap. The \F, m f ) = |1,-1>, |2,1) levels of 
the hyperfine multiplct of Rb split by luq sa 6.8 GHz were 
used. Almost perfect compensation of the density shift 



by spatially varying Zeeman frequency was achieved. The 
transition frequency varied along the sample axis by few 
tens of Hz. Initial inner state with polarization in the 
x — y plane was prepared by a ^ pulse. It was observed 
that the polarization does not remain in the x — y plane 
during free Larmor precession. This was argued to re- 
sult from spatial segregation of atoms with different z 
spin components. However, the confinement potential 
IJ spin dependence was very weak, and the best esti- 
mate of segregation due to the difference in mechanical 
forces experienced by atoms with different spins was at 
least an order of magnitude longer than the time « 0.2 s 
of the z component buildup. 

We argue below that the phenomena of Ref. (I]] are 
explained by coherent evolution of atoms' internal state 
rather than by mechanical segregation in the gas. The 
observed z component profile is readily accounted for by 
interatomic exchange coupling. The transition frequency 
|jj varies along the sample axis, and a short time after 
precession started a gradient of precession angle builds 
up. Now, consider two interacting atoms with slightly 
different polarization due to spatially varying Larmor fre- 
quency. The exchange interaction of these atoms can be 
described [|ll] as precession of each atom's spin around 
the net spin of both atoms. Since both atoms have trans- 
verse polarization, the precession about a net spin (which 
is also transverse) will move the spins out of the x — y 
plane and both of them will acquire a finite z component. 

Exchange effects can be illustrated by a thought ex- 
periment involving a gas of identical atoms with den- 
sity n and spin 1/2 contained in a box. Take the spin 
polarization s to be purely transverse and the same 
for all atoms. For isotropic exchange coupling Ti = 
H J (wos z (r) + ^s(r) • s(r)) d 3 r the polarization s is uni- 
formly precessing, s(t) = in (coswo^ x + smutty)- Now 
consider a test atom passing through the box with spin 
polarization different from that of the other atoms. The 
test atom spin will experience an effective 'magnetic' field 
B = uj i + \s(t) with the exchange part As(i) giving rise 
to Rabi transitions. In a Larmor coordinate frame rotat- 
ing with frequency uiq about the z axis the effective field 
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is just Xn along the gas polarization s, time-independent 
in this frame. Since s is transverse, Rabi transition will 
generate a z component of the test atom polarization, 
even if initially it was in the x — y plane. 

Before accepting this explanation one needs to discuss 
energy conservation. The probabilities to find the test 
atom in the up and down states after coming out of the 
box differ from those in the initial state, since its z spin 
component changes. This means that the test atom en- 
ergy can change by hu>Q. The total energy of the system, 
however, does not change because the spins coupled by 
exchange precess together around the net spin so that the 
total spin is conserved pT| . The change of the net spin 
z component in the box is equal and opposite to the test 
atom spin change, as required by the energy balance. 

Although everything is consistent with energy conser- 
vation, the energy change of Hloq with loq w 6.8GHz much 
higher than other frequencies in the system Q may ap- 
pear counter-intuitive. The temperature T = 600 nK 
corresponds to ksT/h ~ 10 KHz, the trap frequencies 
are (uij_,lj z ) = (230, 7) Hz. However, the characteristic 
exchange frequency is Xn ~ 140 Hz for the typical den- 
sity n = 2 10 13 cm" 3 (see below). This is much higher 
than the transition frequency broadening estimated from 
the precession decay time to be of order of few Hz [Q. 
This makes the exchange induced Rabi transitions at the 
energy Hloq fully coherent, despite that htoo 3> Xn. 

The length corresponding to one Rabi cycle is 



Vt 

- — = 16 am, 
Xn 



(yr = VZT/~r 



(1) 



This is larger than the sample radius r± — 7.3 (im but 
much smaller than the sample length r z — 240 ^m. Since 
movement of an atom by ~ j^cxch is sufficient for rotat- 
ing the spin by £ and moving it out of the x — y plane, 
the exchange coupling is a viable mechanism for spin re- 
orientation in this system. 

The separation of Rb atoms into a gas sample and a 
test particle in the thought experiment is artificial. The 
atoms in Q| share both roles, by inducing precession on 
each other via exchange coupling. One therefore has to 
consider collective dynamics of atom polarization [§,1). 
The Hamiltonian of Rb atoms in a trap has the form 



H 



\i=l,2 3,k=l,2 



:hjhk '■ d 3 r (2) 



H j = ~^ 2 + U j (r), X jk = ^a jk (3) 



where n,- 



tj}j is the density operator. For the states 



used in Ref. |gj the scattering lengths are (an, a 2 2, 012) = 
(100.9, 95.6, 98.2) a with a the Bohr's radius. 

The polarization of internal states is described by 'spin' 
operators with components given by Pauli matrices 



zx{y,z) 



(r) 



2^ 

j,k 



(4) 



and standard commutation algebra 

[<P(r),sV)] =ie aM s^{r)5{r-r') 



(5) 



of spin density operators. 

The system is deep in the cold collision regime, since 
thermal de Broglie wavelength At = h/mvT — 4000 ao is 
much larger than the scattering lengths aj k . We employ 
the forward scattering approximation also known as the 
random phase approximation p9| . The interaction can 
be rewritten in momentum representation as 



hjh k :dr = ipj,pipk,p'^j.,p"^k, P '" (6) 

p+p'=p"+p"' 

(^j,p + ^k,p^j- p ^k-p' + +^j,p + ^kyjpj,-p' + ^k,~p^j 



where p± = p ± q/2. The first term of (g) accounts for 
the forward scattering process, while the second term de- 
scribes exchange scattering. Identifying the operators in 
(0) with the spin density components (JJ) we obtain 



:n\n\ . — 
:h\fi2 + n 2 ni : 



\{n + 2s z f, :n 2 n 2 : = \(h-Ts z f 



In 2 



2(s z Y +4J+J- +4s" 



(7) 

(8) 



In the spin representation the interaction has the form 
1^2 X jk ■■n 3 h k := fn 2 + Ahs z + SX{s z ) 2 + A i2 s 2 (9) 

with u — An + A22 + A12, A = An — A22, and SX = 
A11 + A22 — 2Ai 2 . Spin dynamics is given by 9 t s = j- [s,H], 
where the commutator can be evaluated with the help of 
the relations (JsJ) - After taking the expectation values 
s = (s) we obtain a generalized Bloch equation 

dts + V • j*= fl x s, n = (ujq + Suj)z + 2Ai 2 s (10) 

withj(r) = — ^■(■0jSjfcVr/'fc)+b-c. the spin current. Here 

8uj(r) = \{Ui - U 2 ) + An + 2SXs z (11) 



To make contact with the discussion in Ref. ]l[ we note 
Ui(2)j the occupation probabilities for the 



that hi±s z 



up and down spin. Combined with the form of A and 5X, 
the frequency 5u)(r) can be rewritten as 

fe(r) = i(C/ 1 -f/ 2 )+2(A 11 -A 12 )n 1 -2(A 22 -A 12 )n 2 (12) 

The first term is the Zeeman frequency shift due to the 
trap field inhomogeneity, while the last two terms (identi- 
cal to Eq.(l) of Ref. |Q]) give the density shift. The term 
2Ai 2 s in the expression (|lC]) for fi(r) representing the ef- 
fect of exchange is not considered in Ref. jy. The role 
of this term is subtle. It drops out from the Bloch equa- 
tion (|o|) for s, since s x s = 0. However, since typically 
2Ai 2 |s| ^> |&i>(r)|, this term should be taken into account 
in the Bloch equation for other spin-related quantities, 
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such as the spin current j. The torque s x j (which in 
general is not along the z axis) makes the spin current 
precess so that the x, y and z components intermix. 

The transport equation for the spin current j is derived 
in a similar fashion . We obtain 



d t j + aVs = n(r) X j - 7 j, 



a = \v\ 



(13) 



where the elastic collision rate 7 = 47ra 2 VTn — 20 Hz 
with parameters of Ref. [Q. The term aVs arises in a 
standard way after retaining angular harmonics of the 
lowest order in the transport equation. In Eq.(|l3|) we ig- 
nored the terms such as sV(£/i + U 2 ) and ni\7(U\ — U 2 ), 
since their magnitude is small (see Ref. Q). In this ap- 
proximation, the spin and density dynamics decouple, in 
agreement with the observation m. 

Since the interaction parameters Xjk for Rb coincide 
within 3%, spin is approximately conserved by elastic 
collisions. In this case, the spin current relaxation rate 7 
is the same as for the average particle momentum, and no 
spin relaxation appears in the Bloch equation (10). Elas- 
tic collisions control spin relaxation indirectly, by making 
spin diffusion constant complex |2jj3j. 



The transport equations 
a one dimensional system 



J) , ([13D can be simplified for 
by averaging over sample 



crossection. Large exchange An ~ lu± Jy] leads to fast 
dynamical averaging of spin polarization in each crossec- 
tion with parameters slowly varying along the sample 
length. In averaging Eqs.(M),([L3h we assume Gaussian 



density profile n(p) 



ne 



-p / r ± 



The averaging of the 



terms in Eqs.(pO|),(p^[) quadratic in density and/or spin 
is performed as J ir\p)d 2 p/ J n(p)d 2 p = |n, where n is 
the peak density. After rescaling all coupling constants 



A 



2 \k 



(14) 



and replacing V by one dimensional d x we obtain trans- 
port equations of Leggett-Rice form p0[ 

d t s + d x j — (lo + Slo(x)) z x s (15) 
d t j + ad x s = ( (lo + Slo(x))z + X 12 s) x j 7 j (16) 
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Slo{x) = i(Ui- U 2 ) + - An + SX s 



(17) 



The coupled dynamics of s and j is nonlinear because of 
the exchange precession torque A12S x j in Eq.([T^). 

In the approximation 8lj (x), 7 <C Ai2« p[ one can sim- 
plify transport equations by performing a gradient ex- 
pansion. We first go to the Larmor frame rotating with 
frequency wo, which eliminates luq from Eqs.([l5|),(^6|). 
Next, ignoring the time derivative dtj in Eq.(16) we solve 



it for j in terms of s and d x s, and substitute the result 
in Eq.(|l5|). This gives the Landau-Lifshitz equation ^ij 

d t s - d x (D 1 (s)d x s) = [ho{x) z - D 2 {s)dls] x s (18) 



D 1 (s) 



7 2 + A 2 2 s 



D 2 {s) 



7 2 + A 2 2 s 



(19) 



It is convenient to nondimensionalize Eq. (|T|) . We rescale 
s by |s| max = in, the frequency Slo(x) and scattering rate 
7 by Ai2|s| max , and choose as length unit 



Al2|s| max _ Ai2?T. / 3m A 1/2 s/3 



vi/a 



2k B Tj 



^exch 



14 pm (20) 



Eq. (|l8j) preserves its form, with D 2 — l/(7 2 + s 2 ) and 
D\ = "/D 2 . The dimcnsionless damping 7 is 



4na V'rn 2 a , . 

7 = — = — = -amv T = 47T— ~ 0.31 (21) 

^Xi 2 n fi A T 



since a as 100 ao and the de Broglie wavelength At = 
h/mvT w 4000 a 0. 

The results of numerical simulation of Eq.(jlJ) are 
shown in Fig|l]. The spatial and temporal behavior is 
similar to that in Ref. [Q: The z component builds up 
~ 0.2 s after precession started and then gradually decays 
to zero along with the oscillating transverse component. 
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FIG. 1. Numerical simulation of Eq 
used: constant density n 
Ai2|s| max = 70 Hz. Spatial variation of the transition fre- 
quency 8lu(x) = — £1 cos(2ttx/L) with £1 = 20/n m 6.37Hz 
and sample size L — 10 3 fim. Top: spin distribution evolu- 
tion at < t < 0.2 sec; Bottom: time dependence at x = 0. 
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Spin precession becomes synchronized in different parts 
of the sample (see Fig.||), due to compensation of the 

transition frequency Slo(x) spatial variation by the ex- 
change field Ai2s(x). In our simulation, synchronization 
takes place independently in the domains with s z > 
and s z < 0. Frequency was evaluated as 



/ = ±<J0/dt, e = arg (« a 



ZSq, 



(22) 



During the first 0.2 s of the z component buildup the 
frequency evolves from initial value / = ^-6uj(x) to a 
constant value ±6 Hz in each domain. 

While precession frequencies become synchronized, the 
phase 9 varies within each domain producing spin flux 
between the domains. Spin density s(x) vanishes at 
the domain boundaries x = ijL — ±250 /im (see 
Figs.|l|,H). The number of synchronized domains and 
domain-specific frequency values in general depend on 
the amplitude and characteristic spatial scale of 8w(x). 



In summary, exchange coupling in a trapped gas leads 
to complex collective dynamics of polarization. Polar- 
ized atoms exert torque on the spin current creating a 
z component profile in the presence of spatially varying 
transition frequency. Results of numerical simulation of 
the dynamics of the Rb system jlj are in agreement with 
observations. Surprisingly, the buildup of the z compo- 
nent is accompanied by synchronization of precession fre- 
quencies. In the inhomogeneous state the sample breaks 
into two or more synchronized domains. Spin relaxation 
is caused by spin currents between the domains. 

Synchronized precession should manifest itself in ex- 
periment as transition frequency locking to one value in 
the entire sample, if it is a single domain. Several syn- 
chronized domains formed within the sample will give 
rise to several plateaus in the transition frequency spa- 
tial dependence. Spin density vanishing between different 
domains should be observable by the spatially resolved 
Ramsey fringes technique of Ref. (]| 
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